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ABSTRACT 

Influenza A viruses counteract the ceUular innate immune response at several steps, including blocking RIG I-dependent activa- 
tion of interferon (IFN) transcription, interferon (IFN) -dependent upregulation of IFN-stimulated genes (ISGs), and the activity 
of various ISG products; the multifunctional NSl protein is responsible for most of these activities. To determine the importance 
of other viral genes in the interplay between the virus and the host IFN response, we characterized populations and selected mu- 
tants of wild-type viruses selected by passage through non-IFN-responsive cells. We reasoned that, by allowing replication to 
occur in the absence of the selection pressure exerted by IFN, the virus could mutate at positions that would normaUy be re- 
stricted and could thus find new optimal sequence solutions. Deep sequencing of selected virus populations and individual virus 
mutants indicated that nonsynonymous mutations occurred at many phylogeneticaUy conserved positions in nearly aU virus 
genes. Most individual mutants selected for further characterization induced IFN and ISGs and were unable to counteract the 
effects of exogenous IFN, yet only one contained a mutation in NSl. The relevance of these mutations for the virus phenotype 
was verified by reverse genetics. Of note, several virus mutants expressing intact NSl proteins exhibited alterations in the M1/M2 
proteins and accumulated large amounts of deleted genomic RNAs but nonetheless replicated to high titers. This suggests that 
the overproduction of IFN inducers by these viruses can override NSl -mediated IFN modulation. Altogether, the results suggest 
that influenza viruses replicating in IFN-competent ceUs have tuned their complete genomes to evade the cellular innate immune 
system and that serial replication in non-IFN-responsive ceUs aUows the virus to relax from these constraints and find a new ge- 
nome consensus within its sequence space. 

IMPORTANCE 

In natural virus infections, the production of interferons leads to an antiviral state in ceUs that effectively Umits virus replica- 
tion. The interferon response places considerable selection pressure on viruses, and they have evolved a variety of ways to evade 
it. Although the influenza virus NSl protein is a powerful interferon antagonist, the contributions of other viral genes to inter- 
feron evasion have not been weU characterized. Here, we examined the effects of alleviating the selection pressure exerted by in- 
terferon by serially passaging influenza viruses in cells unable to respond to interferon. Viruses that grew to high titers had mu- 
tations at many normally conserved positions in nearly all genes and were not restricted to the NSl gene. Our results 
demonstrate that influenza viruses have fine-tuned their entire genomes to evade the interferon response, and by removing in- 
terferon-mediated constraints, viruses can mutate at genome positions normally restricted by the interferon response. 



The influenza A viruses are members of the family Orthomyxo- 
viridae and cause annual epidemics and occasional pandemics 
of respiratory disease with considerable public health and eco- 
nomic impact (1, 2). The viral genome comprises eight single- 
stranded RNA molecules of negative polarity, which are encapsi- 
dated into distinct ribonucleoprotein particles (RNPs) containing 
the viral polymerase and multiple nucleoprotein (NP) monomers 
(3-7). Upon recognition of sialic acid-containing cellular recep- 
tors, the virus enters the cell by endocytosis and the viral genome 
reaches the cytoplasm by hemagglutinin-mediated membrane fu- 
sion at the late endosome (8). Transcription and replication of the 
viral RNPs occur in the nucleus and are mediated by the polymer- 
ase, a heterotrimer composed of the PBl, PB2, and PA subunits. 
For transcription, the polymerase uses a cap-snatching mecha- 
nism to generate viral mRNAs, whereas replication leads to the 
production of progeny RNPs similar to those present in the virion 
(3, 4, 6). Finally, progeny RNPs are exported from the nucleus by 



a CRM 1 -dependent pathway and become incorporated in a spe- 
cific manner into virus particles. These progeny virions bud from 
lipid rafts at the plasma membrane in a process that requires the 
specific action of the Ml and M2 proteins (9-13). 
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As a consequence of influenza infection, the cell mounts an 
innate immune response to restrict virus multiplication. The main 
driver of this response is a group of cytokines collectively named 
type I interferons (IFNs) (14). Expression of IFNs requires the 
detection of virus-specific patterns (pathogen-associated molecu- 
lar patterns [PAMPs] ), which are recognized as nonself by cellular 
sensors (pattern recognition receptors [PRRs]) (15), of which 
RIG-I is relevant during the infection of epithelial cells by influ- 
enza virus (16-18). RIG-I becomes activated by recognition of 
double-stranded RNA (dsRNA) containing a 5' -triphosphate (19, 
20) and signals downstream through an interaction with the mi- 
tochondrial protein IPS-l/MAVS/VISA/Cardif The consequence 
of this downstream signaling is the activation of the transcription 
factors IRF3, NF-kB, and ATF-2/c-Jun, which activate IFN tran- 
scription (reviewed in reference 21). Secreted IFNs act in an auto- 
crine and paracrine fashion by activating the IFN receptor. The 
signal is transmitted by the lak/STAT pathway to the nucleus, 
where it activates the transcription of hundreds of IFN-dependent 
genes (ISGs), many of which have antiviral activity (PKR, Mx, 
OAS, etc.) to help the cell to overcome virus infection and restrict 
virus spread to neighboring uninfected cells (14, 22). 

Viruses have evolved a variety of countermeasures to avoid 
cellular innate immunity and maintain reasonable replication ef- 
ficiency. Some viruses can block the activation of PAMP sensors or 
inhibit downstream signaling, thereby diminishing IFN expres- 
sion. Other viruses counteract IFN activity by targeting various 
stages of signaling downstream from IFN receptor, and, finally, 
many viruses directly inhibit the activity of one or several antiviral 
ISG products (reviewed in reference 14). In addition, many lytic 
viruses diminish the level of transcription and/or translation of all 
cellular genes and hence indirectly downregulate the induction of 
IFNs and ISGs. In the case of influenza A viruses, the nonstruc- 
tural NSl protein is the main IFN-counteracting factor. NSl is a 
small multifunctional protein involved in viral protein synthesis, 
viral RNA replication, and virion production (23-27) that also 
modulates cellular posttranscriptional RNA processing and trans- 
port (28-31) (reviewed in reference 32). Influenza A viruses lack- 
ing NSl protein expression cannot multiply in normal, IFN-com- 
petent cells but are viable in cellular systems deficient in innate 
immunity, although viral yields are still reduced compared to 
those of wild-type virus (33, 34). This indicates that the primary 
but not the only role of NSl is to counteract the cellular IFN 
response. This is achieved through a combination of several NSl- 
host cell interactions, some of which are strain specific (35). These 
include (i) inhibiting cellular transcription elongation and post- 
transcriptional RNA processing (28, 36, 37), (ii) blocking RIG-I 
activation (16-18, 38), (iii) interfering with IFN signaling (39, 40), 
and (iv) directly inhibiting specific ISG products, such as PKR and 
RNaseL (41-43). 

In addition to NSl, other influenza virus genes have been re- 
ported to alter the cellular IFN response; the viral polymerase 
inhibits activation of IPS-l/MAVS, possibly as a consequence of 
PB2 localization at the mitochondria (44, 45), while PB1-F2 pro- 
tein appears to modulate the IFN response, although contradic- 
tory results have been reported (46-49). In this report, we have 
undertaken a nonbiased genetic approach to determine the viral 
genes that normally play a role in modulating the cellular IFN 
response. Our results show that the mutations that appear in virus 
populations upon serial passage in non-IFN-responsive cells not 
only are detected in the main IFN response modulator NSl but 



also occur in essentially every viral gene. Moreover, most of the 
virus mutants identified and subsequently characterized induce 
more IFN than wild-type (wt) virus but do not contain alterations 
in NSl. Interestingly, several of these mutants are particularly 
prone to accumulating deletion-containing RNAs derived from 
the genes encoding the polymerase subunits. 

MATERIALS AND METHODS 

Biological materials. In the course of this study, the following influenza A 
strains were used: A/Victoria/3/75 (VIC), A/Puerto Rico/34 (PR8), and 
ANSI (33). Encephalomyocarditis virus (EMCV) was used for the IFN 
bioassay. The MDCK cell line was purchased from the ATCC. The A549 
cell line (50) was obtained from J. A. Melero. The origin of MDCK- V2 and 
MDCK-V5 cell lines has been described (51), as well as the generation of 
A549/pr(IFN-|3).GFP cells (52, 53). A549 cells stably expressing luciferase 
under an interferon-stimulated response element (ISRE) promoter 
[A549/pr(ISRE).Luc] were a kind gift from G. Adolf, Boehringer Ingelheim, 
Austria, and were engineered to express BVDV/NPro [A549/pr(ISRE). 
Luc-BVDV-NP'°] to render them IRF3 deficient and unable to generate 
IFN (54). 

Antibodies used in these procedures included monoclonal antibodies 
against |3-actin (Sigma), phospho-Akt (Ser473; Cell Signaling Technol- 
ogy), and polyclonal antibodies to ISG56 (Santa Cruz), MxA (Santa 
Cruz), phospho-IRF3 (Ser396; Cell Signaling Technology), Akt (pan; Cell 
Signaling Technology), cleaved caspase-3 (Aspl75; Cell Signaling Tech- 
nology), phospho-eIF2a (Ser51; Cell Signaling Technology), and STATl 
(Cell Signaling Technology). The specific anti-NSl and anti-NP rabbit 
antibodies (55, 56), anti-Ml mouse antibodies (12), and anti-NSl rat 
antibodies (24) were described previously. 

Virological techniques. The influenza virus plaque assay was carried 
out on MDCK-V2 cells as described previously (57). Viral plaques were 
stained with crystal violet or by immunocytochemistry using sheep anti- 
sera raised against purified and disrupted X31 (H3N2) virus (anti-X31; 
Diagnostics Scotland) (58). For serial passage, a virus stock generated by 
three consecutive steps of limiting dilution was used as initial inoculum. 
Infections were carried out in sextuplicate at a multiplicity of infection of 
10"' PFU per cell on either MDCK or MDCK-V2 cultures (10^ cells). 
When complete cytopathic effect was observed, the supernatant was col- 
lected, clarified, and used for a further infection cycle. The titers of each 
serial passage virus ranged between 5X10' and 2X10^ PFU per cell. For 
virus purification, virus supernatants were clarified by centrifugation for 
10 min (10,000 X g, 4''C) and sedimented through a 33% sucrose cushion 
in TNE (50 mM Tris HCl-lOO mM NaCI-10 mM EDTA, pH 7.5) for 1 h 
(28,000 X g, 4°C). The virus pellet was resuspended in TNE and sedi- 
mented again on a 50-to-33% sucrose step gradient in TNE. The virus- 
containing interface was collected, diluted in TNE, and sedimented under 
the same conditions (24). 

IFN activity in the culture supernatants was determined by a cyto- 
pathic effect (CPE) reduction bioassay. Culture supernatants from cells 
infected at 5 PFU/cell were harvested after 24 h postinfection (hpi) and 
centrifuged at 1,500 X g for 10 min to eliminate cellular debris. After UV 
treatment to inactivate residual virus, the supernatants was serially diluted 
2-fold and added to A549/BVDV-Npro cell monolayers for 24 h prior to 
infection with EMCV at 0.05 PFU/cell. Monolayers were fixed after 2 to 3 
days and stained with crystal violet to monitor CPE. The number of wells 
protected were converted to IFN bioassay arbitrary units using an IFN-a 
standard. 

Mutant virus screening. For cell sorting, plates of A549/pr(IFN- 
|3).GFP cells were infected at 0.04 FPU/cell. At 8 hpi, cells were 
trypsinized, resuspended in Mg^^- and Ca^^-free phosphate-buffered sa- 
line (PBS), passed through a 0.4-mm-pore-size filter to obtain single-cell 
suspensions, and sorted in an Influx cell sorter (BD Bioscience). Fluores- 
cence-activated cell sorting (FACS) analysis was performed on A549/ 
pr(IFN-P).GFP cells following tiypsinization to obtain single-cell suspen- 
sions and fixation in PBS-1% formaldehyde. Expression of green 
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fluorescent protein (GFP) was examined using a BD FACScan flow cy- 
tometer, and data were analyzed using Flowjo (Treestar). 

Specific mutations were engineered into recombinant virus genomic 
pHH plasmids derived fi-om the VIC strain using the QuikChange site- 
directed mutagenesis kit fi'om Stratagene as recommended by the manu- 
facturer. The mutations were rescued into infectious virus by standard 
techniques (24, 59). 

Protein analyses. Western blotting and immunofluorescence were 
carried out as described previously (60, 61). For immunofluorescence, 
infected A549 cells were washed with PBS, fixed with 1% paraformalde- 
hyde, and permeabilized with 0.5% Triton X-100 for 5 min. The cells were 
blocked with PBS-3% bovine serum albumin (BSA and incubated with 
primary antibody diluted in PBS-1% BSA for 1 h at room temperature. 
After washing with PBS, coverslips were further incubated with goat anti- 
rabbit or goat anti-rat antibodies bound to Alexa 488 or Alexa 594 fluo- 
rochromes. Coverslips were mounted with Prolong reagent and analyzed 
by confocal microscopy using a Leica TSC SP5 microscope. 

For the luciferase assay, A549/pr(ISRE).Luc or A549/pr(ISRE).Luc- 
BVDV-NP" cells were infected at 5 FPU/cell and were treated 7 h later 
with IFN-a (Roferon A; Roche) at 10* units/ml or left untreated. At 13 h 
postinfection, luciferase expression was determined using a luciferase as- 
say system (Promega). 

RNA analyses. The RNA present in purified virions was purified by 
treatment with 0.5% SDS and 200 ixg/ml proteinase K in TNE for 30 min 
at 37°C followed by extraction with phenol-chloroform-isoamyl alcohol- 
hydroxyquinohne and ethanol precipitation (62). Purified RNA was sep- 
arated by electrophoresis on denaturing 4% polyacrylamide- 6 M urea 
gels and silver stained as described previously (62). Reverse transcription- 
PCR (RT-PCR) for PBl, PB2, and PA segments was performed to deter- 
mine the presence of deletion-containing RNAs. To detect the different 
segments, the oligonucleotide 5 ' -GTCACGTCTCATATT AGTAGAAAC 
AAGG-3' was used in combination with either 5' - AGCAAAAGCAGGT 
CAATTATATTCAATATG -3 ' to detect PB2, 5' - CAAAAGCAGGCAAAC 
CATTTGAA-3' to detect PBl, or 5 ' - GCAAAAGCAGGTACTGATTCGA 
GA-3' to detect PA (virus-specific sequences are underlined). The reverse 
transcription reaction was performed for 30 min at 42°C, and then PGR 
was performed for 30 rounds of 94°C for 30 s, 55°C for 30 s, and 68°C for 
2.5 min with a final extension time of 7 min at 68°C. 

Deep sequencing. Deep sequencing of virus populations was per- 
formed using the lUumina HiSeq2000 sequencer with the sequencing kit, 
version 5. After removal of the four nucleotide bar codes used for multi- 
plexing, more than 4 million single-end 76-nucleotide (nt) reads were 
obtained for each sample. Deep sequencing of selected virus mutants was 
performed with the lUumina FiiSeq2000 sequencer using the HiSeq se- 
quencing kit TruSeq v3. In these cases, more than 20 million paired-end 
46-nt reads were obtained. 

Sequence analysis. For alignment of short reads against the viral ge- 
nome, the BWA algorithm (63) was employed, using the sequence of VIC 
as a reference. In all cases BWA alignment parameters were established to 
allow 5 mismatches and no gaps (n, 5; o, 0). SAMtools (64) were used for 
the management of the files generated by BWA. The following strategy 
was used for detecting genomic deletions. Every 46-nt read was split into 
two 23-nt fragments. Each fragment was then independently aligned 
against the VIC genome with BWA software, allowing up to 1 mismatch 
and no gaps. The alignments of two fragments of a single 46-nt read in the 
same strand of the same viral segment at a distance 38 or more nucleotides 
(rndicatrng a gap of 15 or more nucleotides) were scored as positive and 
retrieved. To obtain other junction-compatible read alignments that do 
not lie in the exact middle of a 46-nt read, the entire process was repeated 
with fragments of 20 to 26 nt, 21 to 25 nt, 22 to 24 nt, 24 to 22 nt, 25 to 21 
nt, and 26 to 20 nt. For each sample, all junction-compatible alignment 
events with a gap size of 15 nt or more and that occurred in at least 2 
different reads were considered for subsequent analysis. 
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FIG 1 Diagram of virus serial passage history. MDCK-V2 or MDCK cell cul- 
tures were serially infected at a multiplicity of infection around 0.00 1 PFU/cell 
with A/Victoria/3/75 virus. The initial infecting virus was generated by three 
steps of endpoint dilution to limit the presence of defective virus particles. Six 
replicates were carried out in each cell line, generating passage series VI to V6 
(in MDCK-V2 cells) and Ml to M6 (in MDCK cells). 



RESULTS 

Most of the experimental evidence regarding how influenza virus 
circumvents the IFN response derives from studies of the pheno- 
type of NSl deletion or point mutants and from protein-protein 
interaction studies involving NSl. In contrast, here we performed 
an unbiased genetic screen to identify the viral genes that influence 
the outcome of IFN-mediated innate immune responses to influ- 
enza virus infections. Instead of increasing the IFN selection pres- 
sure (65), we carried out serial passage of A/Victoria/3/75 (VIC) 
influenza virus in non-IFN-responsive cells to release the virus 
from sequence constraints imposed by the IFN response. A dia- 
gram showing the experimental strategy employed is presented in 
Fig. 1 . A wt influenza virus was passaged serially in either MDCK 
or MDCK-V2 cells, which express the V protein of parainfluenza 
virus type 2 (PIV2) and are thus insensitive to IFN-a/p, since V 
targets STATl for degradation in order to block IFN signaling 
(51). Cells were infected at a low multiplicity of infection (MOI) 
(0.001 PFU/cell) to (i) select for efficiently replicating viruses and 
(ii) prevent the accumulation of defective virus (66). Six indepen- 
dent serial passages were carried out in parallel to increase the 
diversity of potential mutants generated. 

Phenotypic and genetic analyses of passaged virus popula- 
tions. After 10 serial passages, the relative efficiency of replication 
of the various passaged stocks was evaluated. MDCK and 
MDCK-V2 cells were infected at a low MOI (0.001 PFU/cell) with 
the parental virus, ANSI virus, or the passage 10 virus series Ml to 
M6 and VI to V6. At various times after infection, the titer of 
infectious virus was assayed on MDCK-V2 cells and the ratios of 
the maximum titers obtained following replication in MDCK-V2 
cells to those in MDCK cells were determined. The parental virus 
used for serial passage produced similar yields in both cell lines, 
while the ANSI virus replicated to approximately 50-fold-higher 
titers in MDCK-V2 cells than in MDCK cells (Fig. 2A). Viruses 
serially passaged in MDCK cells produced about 2-fold-higher 
titers in MDCK-V2 cells than in MDCK cells, but strikingly, vi- 
ruses that had been selected by passage through MDCK-V2 cells 
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FIG 2 Replication efficiency of serially passaged viruses in MDCK and non-IFN-responsive MDCK cells. (A) Replication kinetics of wt and ANSI viruses in 
MDCK-V2 (gray lines) and MDCK cells (black lines). Cell cultures were infected with the initial virus stock (solid lines) or ANSI virus (dotted lines) at an MOI 
of 0.001. At the indicated times after infection, samples were obtained from the culture supernatant, and the viral titers were determined by plaque assay in 
MDCK-V2 cells. (B) The efficiency of replication of the various viruses of the M and V passage series were evaluated by low-multiplicity infections in MDCK-V2 
and MDCK cells after 10 serial passages. Cultures were infected with either the initial virus stock (PO), one of the Ml to M6 or VI to V6 viruses, or ANSI virus 
at a multiplicity of 0.001 PFU/cell. At various times after infection, samples were collected and titers were determined on MDCK-V2 cells. Data are ratios of 
maximal titers of infections in MDCK-V2 to those in MDCK cells, as shown in panel A (averages and standard deviations for 3 replicates). *, P < 0.05. (C) 
Absolute titers of a representative kinetics experiment. (D) The efficiency of replication of some viruses of the M and V passage series was evaluated by using 
low-multiplicity infections in either MDCK-V5 and MDCK cells after 10 serial passages. Cultures were infected with either the initial virus stock (PO) or the Ml, 
M6, VI, V6, or ANSI virus at a multiplicity of 0.001 PFU/cell. At various times after infection, samples were collected and titers were determined on MDCK-V2 
cells. Data are ratios of maximal titers of infections in MDCK- V2 to those in MDCK cells, as shown in panel A (averages and standard deviations for 3 replicates). 
», P<0.05. 



replicated significantly better in MDCK-V2 cells than MDCK cells 
(approximately 10-fold) (Fig. 2B and C), strongly suggesting that 
some of the viruses within this population do indeed replicate 
better than wt virus in non-IFN-responsive cells but are presum- 
ably selected against in IFN-competent cells. To confirm this, 
some of the passage 10 virus populations of the V and M series 
were used to infect either MDCK or MDCK-V5 cells, a different 
cell line that expresses the V protein of parainfluenza virus type 5 
(PIV5), which also blocks IFN signaling (51). Similar results were 
obtained; specifically, the viruses passaged in MDCK-V2 cells 
yielded higher titers in MDCK-V5 cells than in parental MDCK 
cells, whereas viruses passaged in MDCK cells did not (Fig. 2D). 

The phenotype of the passaged virus populations might be due 
to the presence of many mutations affecting several genes or to a 
small number of mutations concentrated in one viral gene, for 
example, NSl. If the latter scenario was the case, then it would be 
expected that deep sequencing of the viral populations would re- 
veal a greater level of variation in one particular viral gene than the 
others. To determine whether this was the case, we carried out 
deep sequencing of purified RNA from passage 10 virions of the V 
or M series. Under the conditions used, an average of 2,342 reads 



per position were obtained, and no overall change in the consen- 
sus sequence was observed. Moreover, no specific position could 
be identified that showed significant sequence variation in the 
VI -V6 passage series compared to the viruses passaged in MDCK 
cells (data not shown) . These results suggested that many different 
mutations across the virus genome may account for the pheno- 
type observed. 

Screening and characterization of relevant virus popula- 
tions. In view of these sequencing results, we next attempted to 
isolate individual mutant viruses from within the population of 
viruses selected by passage through MDCK/V2 cells that activated 
the IFN induction cascade upon infection. In A549 cells engi- 
neered to express GFP under the control of the IFN promoter 
[A549/pr(IFN-p).GFP cells] (52), GFP expression acts as a marker 
for activation of the IFN induction pathway. These cells were in- 
fected with each of the passage 10 VI -V6 virus populations or the 
Ml MDCK-passaged virus and sorted for GFP expression (the Ml 
MDCK-passaged virus could contain spontaneous mutants af- 
fected in IFN counteraction but was included to set a background 
value). This sorting strategy would reveal mutant viruses able to 
induce high levels of IFN but might not identify mutants unable to 
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TABLE 1 Mutations detected in passaged virus populations'^ 
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15.88 




13.33 


9.98 


\j 






226 


T-C 


NO 


0.25 






17.94 




15.00 


11.42 


V 






367 


G-A 


NO 


0.03 


5.95 












LOW V 






1381 


G-A 


NO 


0.07 


16.22 












V 




M 


147 


C-T 


A4lV 


0.41 


8.53 






35.16 


12.19 


35.94* 


V 






309 


G-A 


R95K 


0.33 


46.43- 












V 






370 


A-G 


I115M 


0.27 












6.64 


C 


0/6607 




441 


C-T 


T139I 


0.18 






27.17 








C 


0/6607 




535 


T-C 


NO 


0.25 




9.38 


20.92 




11.78 




LOWV 


83/4734 




696 


G-A 


S224N 


0.14 










11.82 




V 






866 


C-T 


NO 


0.26 






21.10 








V 




N5 


176 


C-T 


NO 


0.00 










9.41 




— 7— 






253 


C-A 


A76E(NS1) 


0.08 












6.53 


c 


1/6226 




459 


A-G 


1145V (NS1) 


0.09 




9.33 




7.68 






V 






540 


G-A 


E172K(NS1) 


1.07 






10.90 


12.68 






c 


3/6226 




553 


A-T 


N176I{NS1) 
M19L(NS2) 


0.20 


9.50 












c 


7/6226 
9/5988 



" The frequency of specific mutations detected in virus populations passaged in MDCK-V2 cells and sorted for 
GFP expression {populations VSl to VS6) is presented, using 5% as the lower threshold. For comparison, the 



frequency of these mutations in a control population passaged in MDCK cells and sorted for GFP expression 
(MS) is also presented. The conservation of these positions in the influenza sequence database is denoted as 
conserved (C; blue}, lov^' variable (LOW V; light blue; 10 < « < 100), or variable (V; n > 100), and the number 
of instances each mutation appears among the total number of sequences screened is indicated. The mutations 
leading to amino acid change are highlighted in green, and the amino acid change is indicated. 



counteract IFN signaling. The infections were carried out at low 
multiplicity (0.04 PFU/cell) to reduce the possibility of dual infec- 
tion of a single cell. Two sorting strategies were used: 1,000 to 
9,000 bulk positive cells from each passaged population were col- 
lected and plated onto cultures of MDCK-V2 cells to rescue the 
selected population of progeny viruses, or alternatively, single 
positive cells were plated onto MDCK-V2 cells to rescue individ- 
ual mutant viruses (see below). 

The genotypes of the bulk virus populations passaged in 
MDCK-V2 cells and selected by cell sorting (populations VSl to 
VS6) were analyzed by deep sequencing of virion RNA on the 
lUumina platform, using viruses from one of the populations pas- 
saged in MDCK cells and sorted in parallel after infection of A549/ 
pr(IFN-3).GFP cells as a control (population MS). The average 
sequence variability per genome position in all populations was 
0.35%. Therefore, we chose 5% (a > 10-fold increase in variabil- 
ity) as the limit to identify potential nucleotide changes in the 
selected population of viruses that may influence the ability of 
influenza A virus to circumvent the IFN response. A number of 



positions in the genome of these selected populations (popula- 
tions VSl to VS6) showed a significant increase in sequence vari- 
ation using this criterion and were not particularly variable in the 
control selected MS population that was passaged in MDCK cells 
(Table 1). The 41 positions that were specifically variable in pop- 
ulations VSl to VS6 mapped to all viral RNA segments except NP 
and identified a potential landscape of mutations associated with 
an enhanced ability to activate the IFN response or a diminished 
ability to inhibit it. Among sequence positions that were specifi- 
cally variable in populations VI to V6, 1 1 showed high variability 
in more than one passage series (Table 1). This suggests that they 
mutated independently during passaging and hence may repre- 
sent positions relevant for the regulation of IFN induction. In 
addition, 24 of the variable positions showed changes that would 
lead to 25 amino acid substitutions in viral proteins (Table 1). 
When the nucleotide changes detected in these 41 variable posi- 
tions were compared to the influenza virus sequence databases, 
many of the synonymous changes corresponded to positions 
that are not phylogenetically conserved (13 out of 17) (Table 1); 
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FIG 3 Replication efficiency of wt and individual mutant viruses in MDCK 
and MDCK-V2 cells. (A) The efficiency of replication of wt and various mu- 
tant viruses was evaluated by determination of maximal titers after low-mul- 
tiplicity infections in either MDCK-V2 and MDCK cells, as indicated in Fig. 2. 



however, most of the amino acid substitutions were observed at 
positions that are hightly phylogenetically conserved (18 out of 
25) (Table 1). These observations suggest that the many synony- 
mous changes imply random drift, whereas the nonsynonymous 
changes were positively selected for during passage through 
MDCK-V2 cells, with some of these mutations resulting in viruses 
that more readily activate the IFN induction cascade, since they 
also induced expression of GFP in the A549/pr(IFN-p).GFP cells. 
As this analysis was done at the population level, it is not possible 
to know whether some of the synonymous changes were cose- 
lected with those involving amino acid alterations. 

Identification and phenotypic characterization of individual 
virus mutants. To isolate individual virus mutants potentially af- 
fected in their activation of the IFN response, single cells scored as 
positive for GFP expression after infection of A549/pr(IFN- 
p).GFP cells with the passage 10 VI to V6 viruses at a low MOI 
(0.04 PFU/cell) were plated individually onto microcultures of 
MDCK-V2 cells, and the supernatants of cells showing a cyto- 
pathic effect were collected. Out of around 1,200 wells containing 
individual GFP-positive cells, viruses were recovered from 109 
wells after endpoint dilution. There are a variety of possible rea- 
sons why not all GFP-expressing cells yield infectious progeny. For 
example, it has been shown that every virus population contains 
many mutants that are not viable but can express virus proteins or 
even replicate intracellularly (67); alternatively, some infected 
cells may be coinfected with defective particles that interfere with 
the replication of wt viruses. Around 75% of these viable viruses 
showed GFP expression higher than that of wt virus when tested 
by reinfection of A549/pr(IFN-3).GFP cells, and 6 of these were 
chosen for further analysis after preliminary phenotype character- 
ization. First, the relative virus yield in IFN-responsive versus 
non-IFN-responsive cells was measured (Fig. 3). All selected mu- 
tants replicated efficiently in multicycle experiments in non-IFN- 
responsive cells. In comparison, in IFN-responsive cells, mutant 2 
grew to much lower titers (Fig. 3A), while mutants 18, 20, and 21 
yielded smaller reductions or slower kinetics (Fig. 3B). Mutants 14 
and 17 gave similar final virus yields in both cell types, although 
mutant 14 showed protracted kinetics in IFN-competent cells 
(data not shown). Comparison of virus plaque size in MDCK 
versus MDCK-V2 cells revealed that only mutant 2 showed a re- 
duction in plaque size in IFN-responsive cells, which was compa- 
rable to that of ANSI virus (Fig. 3C). Furthermore, analyses of 
virus protein expression and localization during high-multiplicity 
infections in A549 cells showed only minor differences from wt 
virus (data not shown). 

Next we analyzed the interplay of these mutants with the in- 
fected cell. The expression of GFP after high-multiplicity infection 
of A549/pr(IFN-P).GFP cells with the 6 selected mutant viruses is 
shown in Fig. 4A. Consistent with these results, all mutant viruses 
triggered the expression of antiviral factors to levels 10- to 50-fold 



Values are averages and standard deviations from 3 experiments. *,P < 0.05; 
**, P < 0.01. (B) Kinetics of virus multiplication at low MOI. The titers at 
various times postinfection of MDCK-V2 (black lines) or MDCK (gray lines) 
cells were determined by plaque assay in MDCK-V2 cells. (C) Plaque size 
analysis of wt or mutant viruses. Wt or mutant viruses (mutant 2 or ANSI) 
were used for plaque assay in parallel in MDCK-V2 or MDCK cell cultures. 
After staining, the plaque size was determined using Image J software. The 
numbers of plaques measured were 176 (wt), 229 (ANSI), and 116 (mutant 2) 
in each cell line. **, P < 0.01; ***, P < 0.001; P < 0.0001. 
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FIG 4 Induction of interferon by wt or mutant influenza viruses. (A) Virus 
progeny from individual A549/pr(IFN-p).GFP cells that were sorted as GFP 
positive and infected at low multiplicity was recovered after limiting-dilution 
infections in MDCK-V2 cells. The ability of these viruses to induce GFP ex- 
pression after infection of A549/pr{IFN-(3).GFP cells was quantified by FAGS, 
using wt and ANSI viruses as references. The results are representative of 3 
independent experiments, in which 10,000 events were measured for each 
sample. (B) The ability of mutant viruses to induce the secretion of antiviral 
factors was tested by infecting of A549 cells at 5 PFU/cell, using wt and ANSI 
viruses as references. Culture supernatants were collected at 24 hpi, and their 
ability to interfere with EMCV infection of A549/BVDV-Npro cells was deter- 
mined by endpoint dilution, as described in Ivlaterials and Methods. •***, p < 
0.0001. 



higher than wt virus after high-multiphcity infection of human 
epithelial A549 cells, although this was lower than the induction 
observed following a parallel infection with ANSI virus (Fig. 4B). 
Likewise, all mutant viruses were able to induce the expression of 
IFN-inducible genes, as demonstrated by the accumulation of lu- 
ciferase reporter upon high-multiplicity infection of A549 ISRE- 
luc cells, in contrast to wt virus infection, which does not activate 
luciferase expression (Fig. 5A). In addition, some of the mutant 
viruses (mutants 2,14, and 2 1 ) had a reduced ability, compared to 
wt virus, to block the upregulation of luciferase by exogenously 
added IFN (Fig. 5B). 

To further analyze the interplay of these mutant viruses with 
the host cell, the expression of ISGs and the induction of specific 
signaling pathways were determined by Western blotting follow- 
ing high -multiplicity infection, using wt and ANSI viruses as ref- 
erences (Fig. 6). These studies revealed the grouping of the mutant 
viruses into three phenotypes. Infection with mutant 2 led to a 
response somewhat similar to that observed in ANSI virus infec- 
tion: strong activation of IRF3, phosphorylation of eIF2a (a 
marker of PKR activity), and partial induction of apoptosis. In 
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FIG 5 Activation of the ISRE by wt and mutant viruses. Cultures of A549/ 
pr(ISRE).Luc (A) or A549/pr(ISRE).Luc-Npro (B) cells were infected at an 
MOI of 5 PFU/cell or mock infected with the indicated mutants or wt virus as 
a reference. A549 NPro/ISRE-luc cells constitutively express the BVDV NPro 
protein, which degrades IRF3; as a result, these cells cannot produce endoge- 
nous IFN (54). At 7 hpi, IFN was added to the indicated samples, and at 13 hpi, 
total cell extracts were prepared and luciferase activity was determined as de- 
scribed in Materials and Methods. Virus infection was verified by determina- 
tion of NP accumulation by Western blotting, using p-actrn as a loading con- 
trol (C). *, P < 0.05; **, P < 0.01; ***, P < 0.001. 



addition, mutant 2 infection showed induction of MxA, which 
was not observed for ANSI virus infection (Fig. 6), probably as a 
consequence of ANSI inducing apoptosis, leading to subsequent 
degradation of STATl by caspases (68). Mutants 20 and 21 be- 
haved similarly; their infection led to partial activation of IRF3, 
expression of ISGs (MxA and ISG56), and activation of Akt (a 
consequence of NSl -dependent activation of PI3K [69]) (Fig. 6), 
but neither of these mutants induced apoptosis significantly. Of 
note, no strict correlation was observed between the levels of ex- 
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FIG 6 Activation of the interferon response by wt or mutant influenza virus 
infection. Tlie induction of ISGs (MxA and ISG56) and activation of cell sig- 
naling pathways (p-IRF3, p-Akt, c-Casp3, STATl, and p-eIF2a) were analyzed 
by Western blotting with specific antibodies at 16 hpi with the indicated mu- 
tant viruses, using uninfected cells (UI) and cells infected with wt or ANSI 
viruses as references. Virus infection was detected by Western blotting with 
antibodies specific for NSl and NP, using p-actin as a loading control. 



pression of ISGs and the corresponding levels of IFN activity se- 
creted (Fig. 4B). These differences may be a consequence of the 
different kinetics of activation of these proteins. Finally, mutants 
14,17, andl8 did not show a clear phenotype in these assays, other 
than robust activation of Akt by infection with mutant 14. 

Genotype of individual virus mutants with an altered re- 
sponse to interferon. To identify the genetic alterations present in 
the mutants described above that could be responsible for their 
observed phenotypes, their full genome sequence was determined 
by deep sequencing of virion RNA. The results are presented in 
Table 2. Nucleotide changes observed in the consensus sequence 
of the mutants compared with the initial virus used to start serial 
passages are in red. In addition, nucleotide changes that, although 
not altering the consensus sequence, reflect a significant propor- 
tion of the virus sequences at the corresponding position (i.e., 
higher than 30%) are in blue. The proportion of the alternative 
nucleotide in all positions recorded (percent alternative nucleo- 
tides) and the genomic heterogeneity found in the original virus 
stock (percent variation in the initial, nonpassaged virus [Mi]) 
and in the virus isolated after sorting of cells infected with the 
control virus populations (percent variation in MS) are also pre- 
sented as a reference. Some of the recorded positions were very 
variable in the control virus populations and are highlighted in 
yellow, as for example mutant 2, position HA 700. We infer that 
these positions in the virus sequence can easily drift and therefore 
consider them probably irrelevant for the phenotype of the corre- 
sponding mutant. Other positions have very low variability in the 
control virus populations (see for instance mutant 2, position NS 
217); these are highlighted in orange and were considered poten- 
tially relevant for the virus phenotype. A number of these nucleo- 
tide changes led to mutations at the protein level, and the amino 
acid changes are indicated (Table 2). Some corresponded to posi- 



tions that had low variation in the control viruses and that were 
phylogeneticaUy conserved in the influenza A database (high- 
lighted in blue) or showed low phylogenetic variation (highlighted 
in light blue). Sequence analysis of virus populations passaged in 
non-IFN-responsive cells and selected by cell sorting showed that 
most of the amino acid changes detected did not map to the NSl 
protein. In fact, the only NSl mutation observed was I64T (NS 
217) in the NSl of mutant 2; none of the other mutant viruses 
encoded altered NSl proteins. The phenotype of mutant 2 was 
similar to that of ANSI virus, for which it is tempting to speculate 
that mutation NSl I64T is responsible. The relevance of the mu- 
tations identified in mutants 14 to 21 is not clear at this point, 
although it is noteworthy that mutants 1 7 to 2 1 contain alterations 
in either Ml (mutant 18, I115M, and mutants 17, 20, and 21, 
D30N) or M2 (mutants 17, 20, and 21, A86S). The mutations 
Ml-Il 15M and M2-A86S have been detected very seldom (either 
0/6,607 or 1/6,456-fold) (Table 2) and could be relevant for the 
mutant phenotype. Mutant 14 has a mutation in the polymerase 
domain of PBl (N455D) (Table 2) that could alter its RNA repli- 
cation efficiency and be the reason for its protracted replication 
kinetics. 

To verify whether the mutations identified in this screen were 
relevant for the observed phenotype, some of them were intro- 
duced into infectious virus by reverse genetics using the VIC virus 
genetic background. Specifically, rescue of mutations NS1-I164T, 
M2- A86S, M 1 -II 1 5M, and M 1 -D30N was attempted. AH rescued 
viruses were viable except mutant M1-D30N, which could be res- 
cued only in combination with mutation M2-A86S (i.e., as found 
in mutants 17, 20, and 21, described above). To analyze the phe- 
notypes of these recombinant viruses, the abilities of rescued vi- 
ruses to activate expression of an ISRE-responsive luciferase re- 
porter and to induce IFN secretion from infected cells were 
examined (Fig. 7). Recombinant mutants NS1-I64T, M1-I115M, 
and M1-D30N/M2-A86S showed luciferase inductions similar to 
those observed in the original mutants 2, 20, and 21, respectively, 
while mutation M2-A86S alone did not show a clear phenotype. 
Similar results were obtained for IFN induction (Fig. 7B). 

Some IFN-inducing virus mutants are particularly prone to 
generating defective interfering particles. Since the induction of 
IFN has been associated with the presence of defective interfering 
(DI) viruses in both influenza viruses and paramyxoviruses (70- 
73), we used bioinformatic analysis of the deep-sequencing data 
for each mutant virus to detect and quantitate the presence of 
internal deletion DI RNAs in the purified virions. An example of 
the results is presented in Fig. 8, the summary of deletion-contain- 
ing RNA accumulations is shown in Tables 3 and 4 and the com- 
plete data set is available via the supplemental material. A large 
number of deletions were observed that were individually very 
infrequent and widely distributed across all genes of the various 
viruses (although they are clearly more abundant in the polymer- 
ase genes) (Table 3). These deletions could be interpreted as the 
result of common polymerase errors during replication. However, 
they would not be efficiently amplified and/or encapsidated, and 
therefore each one was identified only a few times in the data set. 
More abundant deletions were identified that were preferentially 
derived from a subset of viral genes and were differentially repre- 
sented in the wt and mutant viruses (Table 4). Among the abun- 
dant deletions, two size classes could be distinguished: short (100 
to 200 nt) (Table 4; also, see Fig. 8 as an example) or long (almost 
2,000 nt) (Table 4; also, see Fig. 8 as an example). The abundant 



May 2014 Volume 88 Number 9 



jvi.asm.org 4639 



Perez-Cidoncha et al. 



TABLE 2 Mutations detected in individual virus mutants'^ 



MUTANT 


RNA 


POSITION 


NT 


AA 


% 

Alternative 
NT 


% 

Variation 
in Mi 


% 

Variation 
in MS 


CONSERVATION 


2 


HA 


700 


C-A 


T224I 


0.43 


39.88 


47.21 


V 








858 


T-A 


L277IVI 


0.51 


43.1 


44.69 


V 








1296 


A-C 


I423L 


0.43 


12.63 


0.13 


C 


0/3222 




NP 


1302 


A-T 


- 


0.59 


1.03 


0.44 








NA 


1406 


G-A 


D463N 


30.46 


20.12 


3.12 


LOW V 


79/4312 




NS 


217 


T-C 


I64T 


0.59 


0.86 


0.2 


C 


0/6226 


14 


PB2 


1167 


A-G 


- 


9.8 


2.16 


0.3 








PB1 


426 


C-T 


■ 


10.95 


3.5 


0 










1387 


A-G 


N455D 


10.02 


1.01 


0.12 


c 


31/6199 




HA 


132 


C-T 


H35Y 


6.58 


3.24 


0.17 


c 


0/3226 






266 


A-G 


I79M 


8.83 


0.6 


0 


V 


0/3226 




M 


394 


A-C 


- 


17.81 


17.41 


1.04 






17 


PB2 


2057 


AC 


E677A 


42.12 


1.50 


0.08 


c 


3/6100 




HA 


700 


C-A 


T224K 


30.56 


39.88 


47.21 


V 








858 


T-A 


L277M 


35.22 


43.1 


44.69 


V 








1296 


A-C 


I423L 


33.95 


12.63 


0.13 


c 


0/3226 




NP 


891 


A-C 


- 


31.52 


12.25 


24.54 








M 


113 


G-A 


D30N Ml 


43.74 


8.92 


39.63 


V 








394 


A-C 


- 


40.96 


17.41 


1.04 










969 


G-T 


AB6S M2 


44.69 


2.91 


0.14 


c 


1/6454 


18 


HA 


700 


C-A 


T224K 


4.75 


39.88 


47.21 


V 








858 


T-A 


L277M 


4.05 


43.1 


44.69 


V 








1427 


A-T 


E466D 


4.6 


11.61 


17.12 


c 


9/3226 




M 


370 


A-G 


I11SM Ml 


28.49 


0.59 


0.27 


c 


0/6607 






869 


T-C 


- 


29.61 


30.13 


25.95 






20 


HA 


1194 


G-T 


A389S 


32.22 


2.03 


0.19 


c 


1/6326 






1241 


G-A 


- 


30.37 


8.88 


2.57 








NA 


1381 


G-A 




30.94 


2.95 


0.07 


V 








1406 


G-A 


D463N 


38.81 


20.12 


3.12 


LOW V 


79/4312 




M 


113 


G-A 


D30N Ml 


41.71 


8.92 


39.63 


V 








394 


A-C 




45.51 


17.41 


1.04 










869 


T-C 




45.93 


30.13 


25.95 










969 


G-T 


A86S M2 


46.73 


2.91 


0.14 


c 


1/6454 




NS 


176 


C-T 




40.73 


4.39 


0 


V 




21 


M 


113 


G-A 


D30N M1 


38.41 


8.92 


39.63 


V 








394 


A-C 




40.47 


17.41 


1.04 










969 


G-T 


AB6S M2 


36.41 


2.91 


0.14 


c 


1/6454 



" The mutations detected in individual virus mutants are presented. Mutations that appear as such in tlic 
consensus sequence are in red. Those that do not appear as mutations in the consensus sequence but are 
present in a proportion over 30% of the virus sample are indicated in blue. Amino acid changes in yellow are 
very variable in the control virus populations, and those in orange have very low variability in the control virus 
populations. For comparison, the frequencies of these mutations in the initial, nonpassaged virus (Mi) and in a 
control population passaged in MDCK cells and sorted for GFP expression (MS) are also presented. The 
conservation of these positions in the influenza virus sequence database is denoted as conserved (C; blue), low 
variable (LOW V; light blue; 10 < n < 100), or variable (V; n > 100), and the number of instances each 
mutation appears among the total number of sequences screened is indicated. 



RNAs containing short deletions were detected exclusively from 
the PB2 segment of mutants 18 and 20 and were not detectable in 
other segments or mutants (Table 4); the origin of these differ- 
ences is unclear at present. The abundant RNAs containing long 
deletions were exclusively derived from the polymerase genes 
(Table 4 and Fig. 8; also, see the supplemental material) and could 
represent classical DI RNAs (74). The accumulation of long-dele- 
tion RNAs was many- fold higher in mutants 18,20,and21 than in 
the other mutants or the wt virus, as determined by analysis of the 
deep-sequencing data (Table 4; also, see the supplemental mate- 
rial), with mutant 17 representing an intermediate situation. The 
accumulation of deletion-containing RNAs by specific mutants 
was verified by PGR amplification using terminal oligonucleotide 
primers specific for the polymerase genes (Fig. 9). Although prep- 
arations of all mutants were prepared in parallel at low MOT, it 
could be argued that these disparities arose from different passage 
histories of the various mutants. To rule out such a possibility, 
mutants 18, 19, and 20 were plated on MDCK-V2 cells, individual 
plaques were picked, and the viruses were tested for the presence 



of deletion-containing RNAs as described above. The results indi- 
cated that these mutants generated deletion-containing RNAs 
very frequently (Fig. 9, lanes 18.3, 20.1, and 21.7). Interestingly, 
the sizes of the DI RNAs observed after plaque purification were 
different from those that accumulated in the corresponding (non- 
plaque-purified) mutant virus stock, indicating that they arose as 
new deletions in the single plaques and were not carried over from 
the original virus. 

DISCUSSION 

Complex interaction between influenza virus and the ceUular 
innate immune response. Like many other viruses, influenza A 
viruses need to counteract the cellular innate immune response in 
order to achieve a productive infection. Most of this viral coun- 
teraction is mediated by the NSl protein (21, 32, 33, 75), although 
PB2 has also been implicated (44, 45). The approach taken here 
naively questioned the genetic solutions that the virus could adopt 
within its sequence space to optimize its replication under non- 
IFN-responsive conditions and hence could identify other con- 
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FIG 7 Phenotype of recombinant virus mutants. Ttie capacity of mutant vi- 
ruses to induce luciferase under the control of an ISG promoter (A) was tested 
by high- multiplicity infection of A549/pr(ISRE).Luc, using wt virus as refer- 
ence. The capacity of mutant viruses to induce the secretion of antiviral factors 
(B) was tested in infections of A549 cells at 5 PFU/cell, using wt and ANSI 
viruses as references. Culture supernatants were collected at 24 hpi and their 
ability to interfere with EMCV infection of A549/BVDV-Npro cells was deter- 
mined by endpoint dilution and is presented in arbitrary units. *, P < 0.05; 
**»,P< 0.001. 



straints imposed by the innate immune system. The mutations 
identified in several virus populations that replicate to higher 
titers in non-IFN-responsive cells than in IFN-competent cells 
presumably include genetic solutions which are normally re- 
stricted in an IFN-responsive system but which viruses are allowed 
to adopt when the innate immune response is ablated. The exper- 
imental approach used does not necessarily guarantee that iden- 
tified mutations in the selected virus populations affect positions 
normally involved in restricting the IFN response. However, two 
observations suggest that this is likely: (i) the repetitive identifica- 
tion of a fraction of these mutations in several virus populations 
that evolved independently and (ii) the observation that most of 
the amino acid changes detected affect very phylogenetically con- 
served positions. In addition, a number of mutations were iden- 
tified in individual IFN hyperinducer virus clones selected in par- 
allel. The complete set of mutations is summarized in Fig. 10. It is 



Mutant 20 
1 



Copy Number 



2341. 



WT 
1 



2341 



FIG 8 Deletion-containing virus RNAs detected by deep sequencing. The 
figure shows an example of the results obtained by bioinformatics analysis of 
the deep-sequencing data of wt and mutant virion RNAs (PBl RNA segment; 
comparison of mutant 20 and wt virus). The bars indicate the sequence dele- 
tions mapped in the RNA segment coordinates. The gray scale indicates the 
frequency of each particular deletion in the deep-sequencing data. Only dele- 
tions larger than 15 nt and present more than 5 times are shown. The complete 
deletion data set for all viruses is available via the supplemental material. 



clear that amino acid changes were not restricted to the NSl pro- 
tein. Rather, the HA, M 1 , and polymerase genes showed a number 
of mutations comparable to that observed in NSl, while NP and 
NA were the least mutated genes. Among the amino acid changes 
detected in the polymerase, PB2-G693E stands out as potentially 
altering the nuclear localization signal (NLS) and hence importin- 
alpha binding (76), while PA-L589I is located at the PBl binding 
site (77, 78). Two amino acid changes were observed at positions 
phylogenetically conserved within the NSl protein effector do- 
main; E172K and N176I locahze to the cleavage and polyadenyla- 
tion specificity factor (CPSF) -binding region (35, 79-82). Finally, 
M1-T139I affects a conserved position in the protein sequence 
that is located at the dimer interaction surface in the Ml crystal 
structure (83, 84). However, this genomic analysis cannot discern 
whether these mutations affect directly virus counteraction of 
IFN, since they could also reflect random passenger or compen- 
satory mutations that might be linked to the relevant ones. Alto- 
gether, the landscape of mutations suggests that the virus has 
found multiple solutions in its adaptation to efficiently replicate in 
a non-IFN-responsive system, and among the mutants selected by 
IFN-dependent GFP expression, only a fraction had mutations 
that affect the well-known NSl modulator of the innate immune 
response. These results are in contrast with a recently published 
mutagenesis study of the influenza genome (85) in which only the 
hemagglutinin head and NSl protein showed sufficient structural 
flexibility to allow small insertions without altering virus viability. 
Altogether, the results presented here indicate that under normal 
replication conditions, the virus sequence space is restricted in 
essentially every gene by the cellular innate immune response and 
ablating the IFN response allows the virus to explore new possi- 
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TABLE 3 Total accumulation of deletion-containing RNAs in wt and mutant viruses 



No. of total deletion-containing RNAs in segment" 



Virus 


PB2 




PBl 




PA 




HA 




NP 




NA 




M 




NS 




L 


S 


L 


S 


L 


S 


L 


S 


L 


S 


L 


S 


L 


S 


L 


S 


Wt 




2 




5 




2 




2 




27 




8 


4* 


4 


24» 


6 


2 


7 


6 


29 


3 


35 


5 














5* 




35* 




14 


4 


24 


44 


20 


150 


6 








2 








6 


5* 




17 


7 


12 


549 


6 


115 


6 




2 




7 




9 




2 


6* 




18 


756 


37 


120 


254 


831 


16 








2 


2 


4 




21 


2*-l-7 


4 


20 


2,136 


98 


848 


2,698 


3,434 


72 


18 


24 


20 


10 


20 


30 


*6 + 8 


34 


14*-l-4 


18 


21 


676 


15 


564 


18 


1,101 


18 








2 






2* 




5* 





" The number of individual deletion-containing RNAs showing deletions longer than 15 nt that appear at least twice in each sample is given. Data were standardized according to 
the number of deep-sequencing reads obtained for each virus sample and are presented as numbers of short deletions (S; 100 to 200 nt) or long deletions (L; around 2,000 nts). 
Asterisks denote spliced mRNAs from segments M and NS that contaminated the virion preparation. 



bilities compatible with good replication fitness, some of which 
affect levels of IFN induction. 

Poor IFN blockers or IFN hyperinducers? The phenotypic 
and genotypic analyses of the individual virus clones selected 
showed two distinct groups, one represented by clone 2 and the 
other including clones 18, 20, and 21 and to a lesser extent clone 
17. Clone 2 virus is phenotypically similar to ANSI virus, as it 
cannot restrict IFN induction or IFN-induced ISG expression and 
contains an I64T mutation in NSl. This suggested that this single 
change is sufficient to abolish the virus capacity to counteract the 
IFN-mediated response to infection, a contention that was veri- 
fied by the phenotype of the recombinant virus containing NSl- 
I64T as a single mutation (Fig. 7). In contrast, the clone 18/20/21 
group of viruses did not show mutations in NSl protein, and 
hence their capacity to counteract the IFN response should be 
normal. However, these mutants still induced higher levels of IFN 
than wt virus. This apparent contradiction could be explained by 
the fast generation of deleted genomic RNAs in the 18/20/21 
group compared to wt or clone 2 viruses. Whereas all viral clones 
tested generated a low abundance of deletions essentially affecting 
manyviral RNAs, the 18/20/21 group ofviruses accumulated large 
numbers of deletion-containing RNAs originating only from the 
polymerase segments, in line with the early studies characterizing 
influenza DI RNAs (74, 86). However, the clone 18/20/21 group of 
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TABLE 4 Accumulation of abundant deletion-containing RNAs in wt 
and mutant viruses 

No. of abundant deletion-containing RNAs in segment" 



PB2 



PBl 



PA 



HA NP NA M 



NS 



Virus 
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S L 


S 


L S 


LSLSLSLSL 


Wt 












2 










20' 


14 








113 




17 




456 








18 


756 


59 


248 


756 




20 


2,096 


728 


2,632 


3,294 




21 


632 


472 




1,005 





" The number of individual deletion- containing RNAs showing deletions longer than 
15 nt that appear at least 20 times in each sample is given. Data were standardized 
according to the number of deep-sequencing reads obtained for each virus sample and 
are presented as numbers of short deletions (S; 100 to 200 nt) or long deletions (L; 
around 2,000 nt). The asterisk denotes a spliced mRNA from segment NS that 
contaminated the virion preparation. 



PA 
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2323 — 
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726 — 
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Virus genotype 

FIG 9 Experimental verification of the presence of deleted genomic RNAs. 
RNA samples from purified wt or mutant 2, 14, 17, 18, 20, and 21 virions were 
used for RT-PCR using terminal oligonucleotide primers corresponding to the 
PBl, PB2, and PA RNA segments, and the PGR products were analyzed by 
agarose gel electrophoresis. For mutants 18, 20, and 21, individual plaques 
were isolated (plaques 18.3, 20.1, and 21.7), amplified, and used for RT-PCR 
and agarose gel electrophoresis as described above. 
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FIG 10 Mapping the mutations observed in selected virus populations and mutants. The diagram shows the location of mutations detected in virus populations 
and individual virus mutants derived by serial passage in non-IFN-responsive cells and sorting for IFN-dependent GFP expression. The bars show each 
virus-specific protein, and the numbers on the right denote the number of amino acids. Specific features of each protein are indicated. Arrowheads show the 
position of mutations observed in VS virus populations (black) or individual virus mutants (red). Filled arrowheads denote nonsynonymous mutations, whereas 
empty arrowheads indicate synonymous mutations. The arrowhead size indicates the phylogenetic conservation of the mutation site. 



viruses were unable to block IFN-induced ISG expression, and the 
mechanism underlying this defect is still unclear, but it may be a 
result of DI RNA-mediated interference with nondefective virus 
replication. Surprisingly, viral clones 18, 20, and 21 showed mu- 
tations not in the polymerase genes but in the M1/M2 proteins, 
and the relevance of mutations M1-I115M and M1-D30N was 
verified by the phenotype of the corresponding recombinant vi- 



ruses (Fig. 7). One possible explanation is that these mutants are 
particularly prone to the encapsidation of deletion-containing 
RNAs into the progeny virions (11,87). Surprisingly, this does not 
lead to low virus titers (Fig. 3 and data not shown), suggesting that 
deletion-containing RNAs are incorporated in addition to nor- 
mal-segment RNAs. As RIG-I has been shown to preferentially 
associate with short influenza virus RNAs, including DI RNAs 
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from the polymerase segments (70), the properties of 18/20/21 
viruses would suggest that they are particularly good inducers of 
the IFN response due to increased accumulation of deletion-con- 
taining RNAs, rather than poor blockers of IFN induction (since 
they possess intact NSl proteins), in line with the recent report 
showing the importance of DI RNAs in the innate response in vivo 
(88). This phenotype, together with the surprisingly high fitness of 
these viruses, suggests their potential use as attenuated strains. 

In summary, by serial passage of influenza viruses in non- 
IFN-responsive cells, we have obtained virus populations and 
individual virus clones that replicate better in these cells than in 
IFN-responsive cells and are better able to induce IFN and 
IFN-inducible genes than wt virus. The landscape of mutations 
identified indicates that most of the mutant viruses express a wt 
NSl protein and that the mutations are widely distributed 
among all viral genes. This suggests that essentially the entire 
virus gene complement is tuned to counteract the IFN response 
and that it drifts away from such an optimal sequence when the 
IFN response is ablated. 

ACKNOWLEDGMENTS 

We thank G. Adolf, A. Garcia- Sastre, and J. A. Melero for providing bio- 
logical reagents. We thank H. Himmelbauer for advice with deep sequenc- 
ing. The technical assistance of N. Zamarreno and M. Benavides is grate- 
fully acknowledged. 

This work was supported by the Spanish Ministry of Science and In- 
novation (Ministerio de Ciencia e Innovacion) (grant BFU2010-17540/ 
BMC) (J.O.), the Fundacion Marcelino Botin (J.O.), and The Wellcome 
Trust (grant 08775 1/A/08/Z) (R.E.R.). M.P.-C. was a predoctoral fellow 
and awardee of a short-term travel fellowship from the Spanish Research 
Council (CSIC). 

REFERENCES 

1. Neumann G, Noda T, Kawaoka Y. 2009. Emergence and pandemic 
potential of swine-origin HlNl influenza virus. Nature 459:931-939. 
http://dx.doi.org/ 10.1 038/nature08 157. 

2. Molinari NA, Ortega-Sanchez IR, Messonnier ML, Thompson WW, 
Wortley PM, Weintraub E, Bridges CB. 2007. The annual impact of 
seasonal influenza in the US: measuring disease burden and costs. Vaccine 
25:5086-5096. http://dx.doi.Org/10.1016/j.vaccine.2007.03.046. 

3. Neumann G, Brownlee GG, Fodor E, Kawaoka Y. 2004. Orthomyxovi- 
rus replication, transcription, and polyadenylation. Curr. Top. Microbiol. 
Immunol. 283:121-143. http://dxdoi.org/10.1007/978-3-662-06099-5_4. 

4. Resa-Infante P, Jorba N, Coloma R, Ortin J. 2011. The influenza virus 
RNA synthesis machine: advances in its structure and function. RNA Biol. 
8: 1-9. http://dx.doi.Org/10.4161/rna.8. 1 . 1 5302. 

5. Arranz R, Coloma R, Chichon EI, Conesa JJ, Carrascosa JL, Valpuesta 
JM, Ortin J, Martin-Benito J. 2012. The structure of native influenza 
virion ribonucleoproteins. Science 338:1634-1637. http://dx.doi.org/10 
.1126/science.l228172. 

6. Ruigrok RW, Crepin T, Hart DJ, Cusack S. 2010. Towards an atomic 
resolution understanding of the influenza virus replication machineiy. 
Curr. Opin. Struct. Biol. 20:104-1 13. http://dx.doi.Org/10.1016/j.sbi.2009 
.12.007. 

7. Martin-Benito J, Ortin J. 2013. Influenza virus transcription and repli- 
cation. Adv. Virus Res. 87:113-137. http://dx.doi.org/10.1016/B978-0-12 
-407698-3.00004-1. 

8. Palese P, Shaw M. 2007. Orthomyxoviridae: the viruses and their repli- 
cation, p 1647-1689. In Knipe DM, Howley PM (ed), Eields virology, 5th 
ed, vol 1 . Lippincott Williams &. Wilkins, Philadelphia, PA. 

9. Hutchinson EC, von Kirchbach JC, Gog JR, Digard P. 2010. Genome 
packaging in influenza A virus. J. Gen. Virol. 91:313-328. http://dx.doi 
.org/10.1099/vir.0.017608-0. 

10. NodaT, Kawaoka Y. 2010. Structure of influenza virus ribonucleoprotein 
complexes and their packaging into virions. Rev. Med. Virol. 20:380-391. 
http://dx.doi.org/10.1002/rmv.666. 



11. Rossman JS, Lamb RA. 2011. Influenza virus assembly and budding. 
Virology 411:229-236. http://dx.doi.Org/10.1016/j.virol.2010.12.003. 

12. Bourmakina SV, Garcia-Sastre A. 2005. The morphology and composi- 
tion of influenza A virus particles are not affected by low levels of Ml and 
M2 proteins in infected cells. J. Virol. 79:7926-7932. http://dx.doi.org/10 
.1128/JVI.79.12.7926-7932.2005. 

13. Burleigh LM, Calder LJ, Skehel JJ, Steinhauer DA. 2005. Influenza a 
viruses with mutations in the ml helix six domain display a wide variety of 
morphological phenotypes. J. Virol. 79:1262-1270. http://dx.doi.org/10 
.1128/JVI.79.2.1262-1270.2005. 

14. Randall RE, Goodbourn S. 2008. Interferons and viruses: an interplay 
between induction, signalling, antiviral responses and virus countermea- 
sures. J. Gen. Virol. 89:1-47. http://dx.doi.Org/10.1099/vir.0.83391-0. 

15. Takeuchi O, Akira S. 2009. Innate immunity to virus infection. Immunol. 
Rev. 227:75-86. http://dx.doi.Org/10.llll/j.1600-065X.2008.00737.x. 

16. Guo Z, Chen LM, Zeng H, Gomez JA, Plowden J, Fujita T, Katz JM, 
Donis RO, Sambhara S. 2007. NSl protein of influenza A virus inhib- 
its the function of intracytoplasmic pathogen sensor, RIG-I. Am. J. 
Respir. Cefl Mol. Biol. 36:263-269. http://dx.doi.org/10.1165/rcmb 
.2006-0283RC. 

17. Mibayashi M, Martinez-Sobrido L, Loo YM, Cardenas WB, Gale M, Jr, 
Garcia-Sastre A. 2007. Inhibition of retinoic acid-inducible gene I-medi- 
ated induction of beta interferon by the NSl protein of influenza A virus. 
J. Virol. 81:514-524. http://dx.doi.org/10.1128/JVI.01265-06. 

18. Opitz B, Rejaibi A, Dauber B, Eckhard J, Vinzing M, Schmeck B, 
Hippenstiel S, Suttorp N, Wolff T. 2007. IFNbeta induction by influenza 
A virus is mediated by RIG-I which is regulated by the viral NSl protein. 
Cell MicrobioL 9:930-938. http://dx.doi.Org/10.llll/j.1462-5822.2006 
.00841.x. 

19. Cui S, Eisenacher K, Kirchhofer A, Brzozka K, Lammens A, Lammens 
K, Eujita T, Conzelmann KK, Krug A, Hopfner KP. 2008. The C-ter- 
minal regulatory domain is the RNA 5 '-triphosphate sensor of RIG-I. 
Mol. Cell 29:169-179. http://dx.doi.Org/10.1016/j.molcel.2007.10.032. 

20. Takahasi K, Yoneyama M, Nishihori T, Hirai R, Kumeta H, Narita R, 
Gale M, Jr, Inagaki F, Eujita T. 2008. Nonself RNA-sensing mechanism 
of RIG-I helicase and activation of antiviral immune responses. Mol. Cell 
29:428-440. http://dx.doi.Org/10.1016/j.molcel.2007.ll.028. 

21. Wolff T, Ludwig S. 2009. Influenza viruses control the vertebrate type I 
interferon system: factors, mechanisms, and consequences. J. Interferon 
Cytokine Res. 29:549-557. http://dx.doi.org/10.1089/jir.2009.0066. 

22. HaUer O, Staeheh P, Kochs G. 2009. Protective role of interferon- 
induced Mx GTPases against influenza viruses. Rev. Sci. Tech. 28:219- 
231. http://www.oie.int/doc/ged/D6 1 77.PDF. 

23. Falcon AM, Marion RM, Ziircher T, Gomez P, Portela A, Nieto A, 
Ortin J. 2004. Defective RNA replication and late gene expression in 
temperature-sensitive (A/Victoria/3/75) influenza viruses expressing de- 
leted forms of NSl protein. J. Virol. 78:3880-3888. http://dx.doi.org/10 
.1128/JVI.78.8.3880-3888.2004. 

24. Garaigorta U, Falcon AM, Ortin J. 2005. Genetic analysis of influenza 
virus NSl gene: a temperature-sensitive mutant shows defective forma- 
tion of virus particles. J.Virol. 79:15246-15257. http://dx.doi.Org/10.l 128 
/JVI.79.24.15246-15257.2005. 

25. Aragon T, de la Luna S, Novoa I, Carrasco L, Ortin J, Nieto A. 2000. 
Translation factor eIF4GI is a cellular target for NSl protein, a transla- 
tional activator of influenza virus. Mol. Ceil. Biol. 20:6259-6268. http: 
//dx.doi.org/10.1128/MCB.20.17.6259-6268.2000. 

26. Burgui I, Aragon T, Ortin J, Nieto A. 2003. PABPl and eIF4GI associate to 
influenza virus NSl protein in viral mRNA translation initiation complexes. J. 
Gen. Virol. 84:3263-3274. http://dx.doi.Org/10.1099/vir.0.19487-0. 

27. de la Luna S, Fortes P, Beloso A, Ortin J. 1995. Influenza virus NSl 
protein enhances the rate of translation initiation of viral mRNAs. J. Virol. 
69:2427-2433. 

28. Fortes P, Beloso A, Ortin J. 1994. Influenza virus NSl protein inhibits 
pre-mRNA splicing and blocks RNA nucleocytoplasmic transport. EMBO 
J. 13:704-712. 

29. Lu Y, Qian XY, Krug RM. 1994. The influenza virus NSl protein: a novel 
inhibitor of pre-mRNA splicing. Genes Dev. 8:1817-1828. http://dx.doi 
.org/lO.llOl/gad.8.15.1817. 

30. Nemeroff ME, Barabino SM, Li Y, KeUer W, Krug RM. 1998. Influenza 
virus NSl protein interacts with the cellular 30 kDa subunit of CPSF and 
inhibits 3' end formation of cellular pre-mRNAs. Mol. Cell 1:991-1000. 
http://dx.doi.org/10.1016/S1097-2765(00)80099-4. 

3 1 . Satterly N, Tsai PL, van Deursen J, Nussenzveig DR, Wang Y, Faria PA, 



4644 Jvi.asm.org 



Journal of Virology 



Polygenic Nature of Influenza IFN Counteraction 



Levay A, Levy DE, Fontoura BM. 2007. Influenza virus targets the mRNA 
export machinery and the nuclear pore complex. Proc. Natl. Acad. Sci. 
U. S. A. 104:1853-1858. http://dx.doi.org/10.1073/pnas.0610977104. 

32. Hale BG, Randall RE, Ortin J, Jackson D. 2008. The multifunctional NSl 
protein of influenza A viruses. J. Gen. Virol. 89:2359-2376. http://dx.doi 
.org/10.1099/vir.0.2008/004606-0. 

33. Garcia-Sastre A, Egorov A, Matassov D, Brandt S, Levy DE, Durbin JE, 
Palese P, Muster T. 1998. Influenza A virus lacking the NSl gene repli- 
cates in interferon-deficient systems. Virology 252:324-330. http://dx.doi 
.org/10.1006/viro.l998.9508. 

34. Kochs G, Koerner I, Thiel L, Kothlow S, Raspers B, Ruggli N, Sum- 
merfield A, Pavlovic J, Stech J, Staeheli P. 2007. Properties of H7N7 
influenza A virus strain SC35M lacking interferon antagonist NSl in mice 
and chickens. J. Gen. Virol. 88:1403-1409. http://dx.doi.Org/10.1099/vir.0 
.82764-0. 

35. Kochs G, Garcia-Sastre A, Martinez-Sobrido L. 2007. Multiple anti- 
interferon actions of the influenza A virus NSl protein. J. Virol. 81:701 1- 
7021. http://dx.doi.org/10.1128/JVI.02581-06. 

36. Krug RM, Yuan W, Noah DL, Latham AG. 2003. Intracellular warfare 
between human influenza viruses and human cells: the roles of the viral 
NSl protein. Virology 309:181-189. http://dx.doi.org/10.1016/S0042 
-6822(03)00119-3. 

37. Marazzi I, Ho JS, Kim J, Manicassamy B, DeweU S, Albrecht RA, Seibert 
CW, Schaefer U, Jeffrey KL, Prinjha RK, Lee K, Garcia-Sastre A, Roeder 
RG, Tarakhovsky A. 2012. Suppression of the antiviral response by an 
influenza histone mimic. Nature 483:428-433. http://dx.doi.org/10.1038 
/naturel0892. 

38. Gack MU, Albrecht RA, Urano T, Inn KS, Huang IC, Carnero E, Farzan 
M, Inoue S, Jung JU, Garcia-Sastre A. 2009. Influenza A virus NSl targets 
the ubiquitin ligase TRIM25 to evade recognition by the host viral RNA 
sensor RIG-I. Cell Host Microbe 5:439-449. http://dx.doi.0rg/lO.lOl6/j 
.chom.2009.04.006. 

39. Pauli EK, Schmolke M, Wolff T, Viemann D, Roth J, Bode JG, Ludwig 
S. 2008. Influenza A virus inhibits type I IFN signaling via NF-kappaB- 
dependent induction of SOCS-3 expression. PLoS Pathog. 4:el000196. 
http://dx.doi.org/ 10.1371 /journal.ppat. 1 000 1 96. 

40. Pothlichet J, Chignard M, Si-Tahar M. 2008. Cutting edge: innate 
immune response triggered by influenza A virus is negatively regulated 
by SOCSl and SOCS3 through a RIG-I/IFNARl-dependent pathway. J. 
Immunol. 180:2034-2038. http://wv™'.jimmunoLorg/content/180/4/ 
2034.1ong. 

41. Hatada E, Saito S, Fukuda R. 1999. Mutant influenza viruses with a 
defective NSl protein cannot block the activation of PKR in infected cells. 
J. Virol. 73:2425-2433. 

42. Li S, Min JY, Krug RM, Sen GC. 2006. Binding of the influenza A virus 
NSl protein to PKR mediates the inhibition of its activation by either 
PACT or double-stranded RNA. Virology 349:13-21. http://dx.doi.org/10 
.1016/j.virol.2006.01.005. 

43. Min JY, Krug RM. 2006. The primary function of RNA binding by the 
influenza A virus NSl protein in infected cells: inhibiting the 2'-5' oligo 
(A) synthetase/RNase L pathway. Proc. Nafl. Acad. Sci. U. S. A. 103:7100- 
7105. http://dx.doi.org/10.1073/pnas.0602184103. 

44. Graef KM, Vreede FT, Lau YF, McCall AW, Carr SM, Subbarao K, 
Fodor E. 2010. The PB2 subunit of the influenza virus RNA polymerase 
affects virulence by interacting with MAVS and inhibiting IFN-p expres- 
sion. J. Virol. 84:8433-8445. http://dx.doi.org/10.1128/JVI.00879-10. 

45. Iwai A, Shiozaki T, Kawai T, Akira S, Kawaoka Y, Takada A, Kida H, 
Miyazaki T. 2010. Influenza A virus polymerase inhibits type I interferon 
induction by binding to interferon (J promoter stimulator 1 . J. Biol. Chem. 
285:32064-32074. http://dx.doi.org/10.1074/jbc.M110.112458. 

46. Conenello GM, Tisoncik JR, Rosenzweig E, Varga ZT, Palese P, Katze 
MG. 201 1. A single N66S mutation in the PB1-F2 protein of influenza A 
virus increases virulence by inhibiting the early interferon response in 
vivo. J. Virol. 85:652-662. http://dx.doi.org/10.1128/JVI.01987-10. 

47. Dudek SE, Wixler L, Nordhoff C, Nordmann A, Anhlan D, Wbder V, 
Ludwig S. 2011. The influenza virus PB1-F2 protein has interferon- 
antagonistic activity. Biol. Chem. 392:1135-1144. http://dx.doi.org/10 
.1515/BC.2011.174. 

48. Le Goffic R, Bouguyon E, ChevaUer C, Vidic J, Da Costa B, Leymarie O, 
Bourdieu C, Decamps L, Dhorne-PoUet S, Delmas B. 2010. Influenza A 
virus protein PB1-F2 exacerbates IFN-beta expression of human respira- 
tory epithelial ceUs. J. Immunol. 185:4812-4823. http://dx.doi.org/10 
.4049/jimmunol.0903952. 



49. Varga ZT, Ramos I, Hai R, Schmolke M, Garcia-Sastre A, Fernan- 
dez-Sesma A, Palese P. 2011. The influenza virus protein PB1-F2 
inhibits the induction of type i interferon at the level of the MAVS 
adaptor protein. PLoS Pathog. 7:el002067. http://dx.doi.org/10.1371 
/journal.ppat.1002067. 

50. Giard DJ, Aaronson SA, Todaro GJ, Arnstein P, Kersey JH, Dosik H, 
Parks WP. 1973. In vitro cultivation of human tumors: establishment of 
cell lines derived from a series of solid tumors. J. Natl. Cancer Inst. 51: 
1417-1423. 

51. Precious B, Young DF, Andrejeva L, Goodbourn S, Randall RE. 2005. 
In vitro and in vivo specificity of ubiquitination and degradation of STAT 1 
and STAT2 by the V proteins of the paramyxoviruses simian virus 5 and 
human parainfluenza virus type 2. J. Gen. Virol. 86:151-158. http://dx.doi 
.org/10.1099/vir.0.80263-0. 

52. Chen S, Short JA, Young DF, KiUip MJ, Schneider M, Goodbourn S, 
Randall RE. 2010. Heterocellular induction of interferon by negative- 
sense RNA viruses. Virology 407:247-255. http://dx.doi.0rg/lO.lOl6/j 
.virol.2010.08.008. 

53. KUUp MJ, Young DF, Ross CS, Chen S, Goodbourn S, RandaU RE. 201 1. 

Failure to activate the IFN-beta promoter by a paramyxovirus lacking an 
interferon antagonist. Virology 415:39-46. http://dx.doi.0rg/lO.lOl6/j 
.virol.201 1.03.027. 

54. Hilton L, Moganeradj K, Zhang G, Chen YH, Randall RE, McCauley 
JW, Goodbourn S. 2006. The NPro product of bovine viral diarrhea virus 
inhibits DNA binding by interferon regulatory factor 3 and targets it for 
proteasomal degradation. J. Virol. 80:11723-11732. http://dx.doi.org/10 
.1128/JVI.01145-06. 

55. Marion RM, Aragon T, Beloso A, Nieto A, Ortin J. 1997. The N-termi- 
nal half of the influenza virus NSl protein is sufficient for nuclear reten- 
tion of mRNA and enhancement of viral mRNA translation. Nucleic Acids 
Res. 25:4271-4277. http://dx.doi.org/10.1093/nar/25.21.4271. 

56. Coloma R, Valpuesta JM, Arranz R, Carrascosa JL, Ortin J, Martin- 
Benito J. 2009. The structure of a biologically active influenza virus ribo- 
nucleoprotein complex. PLoS Pathog. 5:el000491. http://dx.doi.org/10 
.1371/journal.ppat.l000491. 

57. Ortin J, Martinez C, del Rio L, DavUa M, Lopez-Galindez C, ViUanueva 
N, Domingo E. 1983. Evolution of the nucleotide sequence of influenza 
virus RNA segment 7 during drift of the H3N2 subtype. Gene 23:233-239. 
http://dx.doi.org/10.1016/0378-l 119(83)90055-0. 

58. Hale BG, Kerry PS, Jackson D, Precious BL, Gray A, Killip MJ, Randall 
RE, Russell RJ. 2010. Structural insights into phosphoinositide 3-kinase 
activation by the influenza A virus NSl protein. Proc. Natl. Acad. Sci. 
U. S. A. 107:1954-1959. http://dx.doi.org/10.1073/pnas.0910715107. 

59. Neumann G, Watanabe T, Ito H, Watanabe S, Goto H, Gao P, Hughes 
M, Perez DR, Donis R, Hoffmann E, Hobom G, Kawaoka Y. 1999. 
Generation of influenza A viruses entirely from cloned cDNAs. Proc. Natl. 
Acad. Sci. U. S. A. 96:9345-9350. http://dx.doi.org/10.1073/pnas.96.16 
.9345. 

60. Jorba N, Juarez S, Torreira E, Gastaminza P, Zamarreno N, Albar JP, 
Ortin J. 2008. Analysis of the interaction of influenza virus polymerase 
complex with human cell factors. Proteomics 8:2077-2088. http://dx.doi 
.org/10. 1002/pmic.200700508. 

61. Resa-Infante P, Jorba N, Zamarreno N, Fernandez Y, Juarez S, Ortin J. 
2008. The host-dependent interaction of alpha-importins with influenza 
PB2 polymerase subunit is required for virus RNA replication. PLoS One 
3:e3904. http://dx.doi.org/10.1371/journal.pone.0003904. 

62. de Lucas S, Peredo J, Marion RM, Sanchez C, Ortin J. 2010. Human 
Staufenl protein interacts with influenza virus ribonucleoproteins and is 
required for efficient virus multiplication. J. Virol. 84:7603-7612. http: 
//dx.doi.org/10.1128/JVI.00504-10. 

63. Li H, Durbin R. 2009. Fast and accurate short read alignment with Bur- 
rows-Wheeler transform. Bioinformatics 25:1754—1760. http://dx.doi 
.org/ 1 0. 1 093/bioinformatics/btp324. 

64. Li H, Handsaker B, Wysoker A, FenneU T, Ruan J, Homer N, Marth 
G, Abecasis G, Durbin R. 2009. The sequence alignment/map format 
and SAMtools. Bioinformatics 25:2078-2079. http://dx.doi.org/10 
. 1093/bioinformatics/btp352. 

65. van Wielink R, Harmsen MM, Martens DE, Peelers BP, Wijffels RH, 
Moormann RJ. 2012. Mutations in the M-gene segment can substantially 
increase replication efficiency of NSl deletion influenza A virus in MDCK 
cells. J. Virol. 86:12341-12350. http://dx.doi.org/10.1128/JVI.01725-12. 

66. von Magnus P. 1951. Propagation of the PR8 strain of influenza A virus in 
chick embryos. II. The formation of incomplete virus following inocula- 



May 2014 Volume 88 Number 9 



jvi.asm.org 4645 



Perez-Cidoncha et al. 



tion of large doses of seed virus. Acta Pathol. Microbiol. Scand. 28:278- 
293. 

67. Brooke CB, Ince WL, Wrammert J, Ahmed R, Wilson PC, Bennink JR, 
Yewdell JW. 2013. Most influenza A virions fail to express at least one 
essential viral protein. J. Virol. 87:3155-3162. http://dx.doi.org/10.1128 
/JVI.02284-12. 

68. King P, Goodbourn S. 1998. STATl is inactivated by a caspase. 1. Biol. 
Chem. 273:8699 - 8704. http://dx.doi.org/ 10. 1074/jbc.273. 15.8699. 

69. Hale BG, Jackson D, Chen YH, Lamb RA, Randall RE. 2006. Influenza 
A virus NSl protein binds p85beta and activates phosphatidylinositol-3- 
kinase signaling. Proc. Natl. Acad. Sci. U. S. A. 103:14194-14199. http: 
//dx.doi.org/ 1 0. 1 073/pnas.0606 109103. 

70. Baum A, Sachidanandam R, Garcia-Sastre A. 2010. Preference of RIG-I 
for short viral RNA molecules in infected cells revealed by next-generation 
sequencing. Proc. Natl. Acad. Sci. U. S. A. 107:16303-16308. http://dx.doi 
.org/10.1073/pnas.l005077107. 

71. Strahle L, Garcin D, Kolakofsky D. 2006. Sendai virus defective- 
interfering genomes and the activation of interferon-beta. Virology 351: 
101-111. http://dx.doi.Org/10.1016/j.virol.2006.03.022. 

72. Easton AJ, Scott PD, Edworthy NL, Meng B, Marriott AC, Dimmock 
NJ. 2011. A novel broad-spectrum treatment for respiratory virus infec- 
tions: influenza-based defective interfering virus provides protection 
against pneumovirus infection in vivo. Vaccine 29:2777-2784. http://dx 
.doi.org/10.1016/j.vaccine.2011.01.102. 

73. Killip MJ, Young DF, Gatherer D, Ross CS, Short JA, Davison AJ, 
Goodbourn S, Randall RE. 2013. Deep sequencing analysis of defective 
genomes of parainfluenza virus 5 and their role in interferon induction. J. 
Virol. 87:4798-4807. http://dx.doi.org/10.1128/JVl.03383-12. 

74. Jennings PA, Finch JT, Winter G, Robertson JS. 1983. Does the higher 
order structure of the influenza virus ribonucleoprotein guide sequence 
rearrangements in influenza viral RNA? Cell 34:619-627. http://dx.doi 
.org/10.1016/0092-8674{83)90394-X. 

75. Ehrhardt C, Seyer R, Hrincius ER, Eierhoff T, Wolff T, Ludwig S. 
2010. Interplay between influenza A virus and the innate immune sig- 
naling. Microbes Infect. 12:81-87. http://dx.doi.0rg/lO.lOl6/j.micinf 
.2009.09.007. 

76. Tarendeau F, Boudet J, GuilUgay D, Mas PJ, Bougault CM, Boulo S, 
Baudin F, Ruigrok RW, Daigle N, EUenberg J, Cusack S, Simorre JP, 
Hart DJ. 2007. Structure and nuclear import function of the C-terminal 
domain of influenza virus polymerase PB2 subunit. Nat. Struct. Mol. Biol. 
14:229-233. http://dx.doi.org/10.1038/nsmbl212. 

77. He X, Zhou J, Bartlam M, Zhang R, Ma J, Lou Z, Li X, Li J, Joachimiak 
A, Zeng Z, Ge R, Rao Z, Liu Y. 2008. Crystal structure of the polymerase 



PA(C)-PB1(N) complex from an avian influenza H5N1 virus. Nature 454: 
1123-1126. http://dx.doi.org/10.1038/nature07120. 

78. Obayashi E, Yoshida H, Kawai F, Shibayama N, Kawaguchi A, Nagata 
K, Tame JR, Park SY. 2008. The structural basis for an essential subunit 
interaction in influenza virus RNA polymerase. Nature 454:1127-1131. 
http://dx.doi.org/10.1038/nature07225. 

79. Li Y, Chen ZY, Wang W, Baker CC, Kmg RM. 2001. The 3'-end-processing 
factor CPSF is reqiured for the splicing of single-intron pre-mRNAs in vivo. RNA 
7:920-931. http://dx.doi.org/10.1017/S1355838201010226. 

80. Noah DL, Twu KY, Krug RM. 2003. Cellular antiviral responses against 
influenza A virus are countered at the posttranscriptional level by the viral 
NSIA protein via its binding to a cellular protein required for the 3' end 
processing of cellular pre-mRNAs. Virology 307:386-395. http://dx.doi 
.org/10.1016/S0042-6822(02)00127-7. 

81. Shimizu K, Iguchi A, Gomyou R, Ono Y. 1999. Influenza virus inhibits 
cleavage of the HSP70 pre-mRNAs at the polyadenylation site. Virology 
254:213-219. http://dx.doi.org/10.1006/viro.1998.9555. 

82. Twu KY, Noah DL, Rao P, Kuo RL, Krug RM. 2006. The CPSF30 
binding site on the NSIA protein of influenza A virus is a potential anti- 
viral target. J. Virol. 80:3957-3965. http://dx.doi.Org/10.1128/JVI.80.8 
.3957-3965.2006. 

83. Arzt S, Baudin F, Barge A, Timmins P, Burmeister WP, Ruigrok RW. 
2001. Combined results from solution studies on intact influenza virus Ml 
protein and from a new crystal form of its N-terminal domain show that 
Ml is an elongated monomer. Virology 279:439-446. http://dx.doi.org 
/10.1006/viro.2000.0727. 

84. Sha B, Luo M. 1997. Structure of a bifunctional membrane-RNA binding 
protein, influenza virus matrix protein Ml. Nat. Struct. Biol. 4:239-244. 
http://dx.doi.org/10.1038/nsb0397-239. 

85. Heaton NS, Sachs D, Chen CJ, Hai R, Palese P. 2013. Genome-wide 
mutagenesis of influenza virus reveals unique plasticity of the hemagglu- 
tinin and NSl proteins. Proc. Nad. Acad. Sci. U. S. A. 110:20248-20253. 
http://dx.doi.org/10.1073/pnas.1320524110. 

86. Nayak DP, Chambers TM, Akkina RK. 1985. Defective-interfering (DI) 
RNAs of influenza viruses: origin, structure, expression, and interference. 
Curr. Top. Microbiol. Immunol. 114:103-151. 

87. Rossman JS, Jing X, Leser GP, Balannik V, Pinto LH, Lamb RA. 2010. 
Influenza virus m2 ion channel protein is necessary for filamentous virion 
formation. J. Virol. 84:5078-5088. http://dx.doi.org/10.1128/JVI.00119-10. 

88. Tapia K, Kim WK, Sun Y, Mercado-Lopez X, Dunay E, Wise M, Adu M, 
Lopez CB. 2013. Defective viral genomes arising in vivo provide critical 
danger signals for the triggering of lung antiviral immunity. PLoS Pathog. 
9:el003703. http://dx.doi.org/10.1371/journal.ppat.1003703. 



4646 Jvi.asm.org 



Journal of Virology 



